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s: 

Abstract. We study the correlation function of middle spins, i. e. of spins on 
2 ' intermediate sites between two adjacent parallel lattice axes, of the spatially anisotropic 

t^ . Ising antiferromagnet on the kagome lattice. It is given rigorously by a Toeplitz 

determinant. The large-distance behaviour of this correlation function is obtained 
by analytic methods. For shorter distances we evaluate the Toeplitz determinant 
numerically. The correlation function is found to vanish exactly on a line Jd{T) in 
' ^ I the T ^ J (temperature vs. coupling constant) phase diagram. This disorder line 

divides the phase diagram into two regions. For J < Jd{T) the correlations display the 
QJ , features of an unfrustratcd two-dimensional Ising magnet, whereas for J > Jd{T) the 

correlations between the middle spins are seen to be strongly influenced by the short- 

^S| ■ range antifcrromagnctic order that prevails among the spins of the adjacent lattice 

K*" ' axes. While for J < Jd{T) there is a region with ferrimagnetic long-range order, the 

r^ ' model remains disordered for J > Jd{T) down to T = 0. 
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Figure 1. Anisotropic kagome model (left) and auxiliary model (right). The factors 
of the transfer matrix acting in and between the different rows are indicated on the 
right of the auxiliary model. 



1. Introduction 

Recent experimental work on the mineral volborthite has stimulated several theoretical 
investigations of the Heisenberg antiferromgnet on the spatially anisotropic kagome 
lattice, figured], [I1I21E]- In this model, the exchange coupling in one of the three lattice 
directions (J) differs from the couplings in the two other directions (J'). 

The spin ensemble is thus divided into a set of parallel chains of spins with coupling 
J and a set of middle spins, which are coupled to the chain spins by the exchange 
constant J'. The results for the ground state of this model, obtained in the studies 
[I1E113], differ substantially from each other, particularly in the quasi-one-dimensional 
limit J ^ J' where the spin chains are only weakly coupled to the middle spins . 
In this article, motivated by the difficulties encountered in approximate treatments of 
the Heisenberg model, we address the Ising model with anisotropic couplings on the 
kagome lattice (AIK). As was first demonstrated by Kano and Naya [1], this model can 
be solved exactly. Its thermodynamic properties are discussed in detail in recent work 
by Wei Li et al. [5j. For J ^ J', the ground states are immediately clear. In the region 
J < J', the ground state is ferrimagnetically ordered, i. e. the spins on the chains are 
ordered ferromagnetically, and so are the middle spins, but their direction is antiparallel 
to the chain spins. For J > J', the chain spins are ideally ordered antiferromagnetically 
in the ground state leaving the middle spins completely uncoupled. As a consequence, 
the chains are also decoupled from one another, so that the spins of each chain can form 
either one of the two antiferromagnetic orders. A simple counting argument yields a 
rest entropy of S^^f^^ = |ln2 ^ 0.231 per spin in this case. For the isotropic model, 
J = J', the value of the rest entropy does not follow from a simple counting argument. 
Kano and Naya [1] find in this case 
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S-^f = — — / / In [21 - 4(cos(a;i) + cos(a;2) + cos(a;i + U2))]dujidu2 
247r Jq Jq 

- 0.5018 . (1) 

In the present paper, we shall reformulate the transfer matrix formalism of Kano 
and Naya with the aim to obtain not only the partition function, but also the correlation 
function xi'f^i''^) between two middle spins in the same row (see figure [1]). We expect 
this to be the most interesting of the various spin-spin correlation functions of the AIK 
model. 

The article is organised as follows. In section 2, we consider a slightly modified 
model that has been devised by Kano and Naya to allow for the application of the 
transfer matrix method, and that reproduces the AIK model in a certain limit. We 
present the transfer matrix, discuss its symmetry properties and convert it into a 
product of Grassmann functionals [6]. In section 3, we discuss the structure of the 
eigenvectors of the transfer matrix and calculate its eigenvalues. The latter are then 
used to construct the phase diagram of the AIK model as far as this is possible without 
knowledge of the correlation function. Section 4 is the central part of this article; it 
is devoted to our study of the above mentioned two-spin correlation. As in numerous 
previous investigations of correlations in Ising models on various two-dimensional lattices 
[3, El El [10, [m [121 [131 [HI [Ej, the correlation function considered here is found to 
be represented rigorously by a Toeplitz determinant. By evaluating this determinant 
[I6l[17] we obtain analytic expressions for xij^^ '^) iii the limit of large distances between 
the spins, \n — m\ 3> 1, in various regions and on various special lines of the phase 
diagram. We support our asymptotic results by exact numerical evaluations of the 
Toeplitz determinant, which yield the correlation function also for short distances. In 
section 5 we present and discuss our results. Technicalities of the calculations that lead 
to these results are contained in appendices A, B, C In appendix D, we present a simple 
perturbative calculation, which is valid in the limit of large chain coupling, J / J' ^ 1, 
and which qualitatively reproduces the behaviour of the correlation function in this 
limit. 

After completion of this work, we became aware of a publication by M. B. Geilikman [TS] . 
which addresses the same topic as the present paper. Using the methods of Vdovichenko 
[T9l[20], Geilikman arrives at the same Toeplitz determinant as we, to represent the spin- 
spin correlation function. However, his evaluation of this determinant yields results that 
differ qualitatively from our expressions for xi^^i^) iii the most interesting part of the 
phase diagram. 

2. Transfer matrices 

Following Kano and Naya [1] , we split each intermediate spin of the AIK into two Ising 
variables aj, ctj+i, which are coupled ferromagnetically with a coupling /. At some point 
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of our calculation we take the limit / — ?■ oo so that these two Ising spins freeze into 
the original intermediate spin. The auxiliary model, which results from the splitting, 
is shown in the right part of figure [H Obviously, it is periodic with period 4 in the 
(vertical) transfer direction. We consider a lattice of M (M divisible by 4) rows with 
A'" (N even) sites in each row and impose periodic boundary conditions in both lattice 
directions. For the auxiliary model, the transfer matrix, whose highest eigenvalue yields 
the free energy of this model, is given by (see figure [1]) 



where 



with 



To T2 To Ti To T3 To T,, (2) 



^^^^o<, (3) 

(4) 
(5) 
(6) 



K = J/T , K' = J'/T and tanh(a) = e"^^ . (7) 





N-l 


To 


= (g2i^' _ g-2i^')^/2 11^^ 




n=0 


Ti 


= Q-^ 2^n = <^n ^n + 1 


T2 


= g^Ei^o"^'(-in + l-in + 2-l), 


T3 


= g^Ei^0''/'{-in-2Wl-l), 



(Boltzmann's constant ks is set equal to unity throughout this paper.) 
Instead of K and K' we shall use 

^ = tanh(i^) and $ = tanh(Ji^') (8) 

as variables in the sequel. 

The expression ([2]) is only one of various possibilities to represent the transfer 
matrix of the auxiliary model. In fact, it is an inconvenient representation, since it is 
not symmetric. Instead we can work with the following matrix 

T^^^ To Ti To Tg/^ ■ Tg/^ To Ti Tq T^^^, (9) 

which in turn is equivalent to 

T = T2 To Ti To T3 ■ T3 To Ti To T2 (10) 

in the limit / — t- 00, since T2 and T3 become projectors in this limit, T^ = Ti, i = 
2, 3. T is real and symmetric since T3 = s T2 s, where s shifts all spin variables by one 
lattice unit, s a^ s = 0-^+1^ ^ ^n ^ — ^n+i' ^"^^ therefore commutes with To and T^. 

In the sequel, we follow the transfer-matrix method of Schulz, Mattis and Lieb [9]. 
The necessary calculations are straight-forward. However, since the kagome lattice and 
the auxiliary lattice are non-Bravais lattices, some pecularities arise in comparison to 
the square lattice case. Therefore, we sketch the transfer-matrix procedure briefly here. 
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In order to exhibit the free fermion nature of the AIK model, one transforms the 
expressions in the exponents of Tq, ... T^ into quadratic forms in fermion operators 
bn, h\^ by applying the Jordan- Wigner transform [21] : 

< = 2hih^ - 1, 

hi + 6n and af = {h^ 



<Jo = Uq^-Uo — "n 



K) e-^^ 



^^.7 = J 3 



<n<N. 



This yields 



To 
Ti = e 

T, 



4$ 



N/2 



-Na p(2a+J7r)E„=o f^f-' 



1_$2 



or, after Fourier transformation to momentum space 



(11) 
(12) 



(13) 

(14) 
(15) 
(16) 



To = 
Ti =e- 

^2,3 = e" 



16$^ 



:i-$2)2 

-Khi 



N/A 



,(2a + ivr) Iiq 



-2iN/A „2^^2,3 



with 



/lo 

/i2,3 



^fc ) 



T<A;<7r 

5^ [z sin(A:) {h{hl, + 6,6_,) + cos(A;) {h\h, - h,hl 

7T<.k<.7T 

7r<fc<7r 

- cos(A;) f 6i6, - b.bl ± 6,6,_, ^ 616^^ 



— 7r<A;<7r 



(17) 



(19) 



(20) 



We work with periodic boundary conditions. For states with even (odd) numbers 
of fermions this requires that 6)^ = 6j (fejy = — fep) (cf. reference [5]). To satisfy these 
conditions one chooses in b^ = -4= J2k ^^^^ ^k wave numbers k = ^m or k = ^(m+1/2), 
m integer, for the even and odd states, respectively. 

The non-Bravais structure of the auxiliary lattice in figure [1] manifests itself in the 
structure of the hamiltonians /i2,3: they couple not only fermionic states with wave 
numbers k and —k as in the square lattice case, but also states with wave numbers k 
and ±k + n. To deal with this situation we divide the Brillouin zone — vr < k < tt into 
4 subintervalls, -it < k < -it/2, -7t/2 < k < 0; < k < 7t/2, 7t/2 < k < tt and 
choose iV to be a multiple of 4. Then, /iq, hi and /i2,3 are diagonal in q and can be 
represented in the reduced Brillouin zone < g < tt /2. For instance. 
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ho 


= y, ho, + Kr with 

0<q<n/2 




hoq 


= &l.+g&-.+, + bl^b_^ + b% + bi_^b^_g , 


, and 


i^even 


= bl^b^, + 6l,/2&-V2 + blK + ^l/2^V2 • 





(21) 

The extra term h^^^"' corresponds to the values k = — vr, — vr/2, 0, 7r/2 of the original 
Brillouin zone. It occurs only if ho acts in the space of even states. Then, 

/ 16$^ V^^ T-r 

To-- =7^3^^ To, n ^0. (22) 

where To, = e'^'^'^ol"" and Tog = e'^"'^"'' . Similarly, Ti and T2^3 are products of matrices 
acting in the subspaces of the individual wave numbers q. We remark that while these 
products consist of N/A — 1 factors corresponding to the wavenumbers < g < 7r/2, 
there is only 1 factor of the type Tq^ in each of the matrices Tq, Ti, T2, and T3. In 
the present paper we are always interested in the thermodynamic limit A^ — )■ 00 of our 
results. In this limit the factors Tq, play no role and will therefore be neglected in the 
sequel. 

3. Eigenvalues and eigenvectors of the symmetric transfer matrix 

Our next aim is to determine the largest eigenvalue and the corresponding eigenfunction 
of the symmetric transfer matrix (fTO|) . It consists of factors that are all exponentials 
of quadratic forms of fermion operators. Hence T is also an exponential of a quadratic 
form of fermion operators. To determine this last exponential, the multiplication of 
the individual factors in equation ^U\ has to be performed explicitly. To achieve this, 
we find it convenient to work with the representation of the fermion operators bg, 6j 
by Grassmann variables 7^, 7* In appendix A we briefly describe the steps that are 
necessary to convert operators such as To,l,2,3 into the so called matrix form. For 
details, we refer the reader to the book by F. A. Berezin [22], in particular to chapter 
6. We introduce the shorthand notation 7^ = {7g,7_g,7_^+g,7^_J = {7i, 72 > 73^74} • 
(There is no need to distinguish the row vector 7 from its adjoint column vector 7 * 
here, since dyadic products 7^ ■ 7^* will not occur.) In matrix form, Tq, Ti and T2^3 read 
(here and in the sequel we use the shorthand notation cos(g) = Cq, sin(g) = Sg ) 

To(7*,7)=fT^Ve-^°^.*^., ho = I (23) 



$2/ - ' ■-" $' 



d + ^2)2 _ 4 2^2 ^* h,^ + -7* h[ 7* 7 h'* 7 

Ti(7*, 7) = ^ ' ,,,7 e^' ' ^' 2^' ' ^' 2^" ^' ^-^ (24) 



* 
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with 

- _ (tT^^ \ h' - ^''^^ hi""' M (2^^ 

T,(7*, 7) = \ e^^* ^' + ^'^ ^^ + r^* ^'^ ^^* - r^ ^'^* ^^ (26) 

with 

fc.=f '^ . '''"'V fci=-f /'. "'"If;" \ I (27) 

Here the hmit i ^ oo has already been taken. 
From the remarks after equation (fTO!) of section 2 it follows that 



T = Ti/2 ■ Ti/2, where 

Ti/2 = T2 To Ti To s T2. (28) 

Thus it suffices to determine the eigenvalues and the eigenstates of Ti/2- Following 
the technique described in appendix A, we determine the matrix form of T1/2: 



Ti/2 9(7*7)= / d\*dX djj*djj diy*di/ dp*dp e 



^ e^^ ^ e~^ ^ e~P P 
T2(7*,A) to(A*,/i) Ti(/i*,i/) to(z/*,p) t2(sp*,7) 
^^jli^7. + ^7li/'7l47,i/'*7, (29) 

TV2 =\{T^, (30) 

<? 

(We have suppressed the momentum label q in the first two lines of this equation.) 
The shift operator acts on p* as s p* = {pl, P2, —p%, —pl}- 

We point out here that being quadratic in fermion variables, the exponent in T1/2 g, 
equation ( l29l) . is the hamiltonian of a system of noninteracting fermions. 

The general structure of the hamiltonians H and H' derives from the structure of 
the hamiltonians /lo, hi, h[, /12 and /12, equations ( 123|) . (!25|) . (!27|) : /lo, /^i, /i2 and /I'j^, /i2 
are of the form 

DJl = \ . ^ I and VJl-u, u = \ j , respectively, (31) 

where a, 6, c, and c? are real scalars. Products of Wl are again of the form Wl and 
furthermore [OJt, m] = 0, m^ = 1. Hence, in equation (!29|) . _ff is of the form DJl and iJ' 
is of the form VJt u. Furthermore, since T1/2 is a real symmetric quadratic form in the 
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8 



space of the Grassmann vectors 7 ? H must be symmetric and H' must be antisymmetric. 
Hence H and H' must have the following structure: 



H 



a iba^ 
-iha^ d 



and H' 



a' iV o^ 
d 



(T« 



ib'a' 



-a'' 



.(32) 



The real scalars a, b, d and a', b', d' are rational functions of ^ and $. Explicit expres- 
sions for these quantities are given in appendix B. Their calculation is straight-forward. 

Turning to the eigenvalue problem for the transfer matrix we first consider the space 
of even states. As an even eigenstate we choose in Grassmann representation 



1 



'Pgi.l*) = </'o - i<Pi 7i72 + i<t^2 llll - h -^illll + 7273) 



V2' 



+ 04 7i 72 73 74 • 



(33) 



Again the momentum variable q is suppressed here. Only products of Grassmann 
variables with total momentum g = mod vr are included in this state. The coefficients 
00,1,2,3,4 of the eigenstate ( l33l) are calculated from 



dX*dX ri/2(7*,A) e~^^ 03(A*) = V^^r 0.(7*) 



<g) 



where (see (132 



(34) 



7^1/2(7*, A) =vexp {a(7i Ai + -f2^2) + ^(73^ + llK) 

+ iKll^S - 72^4 - 73^1 + 1*4^2) 

- ^a 7i72 + td 7374 - b (7^74 + 7273) 

— m'Aj^Ag + id'X^X^ + b'^X^X^ + AgAg)}, 



(35) 



and we observe that the products with nonzero (mod tt) momenta, 7i73 and 7274, 
are indeed not generated by the action of T1/2 on 0^(7*). 
Equation flMl) results in 



T 



h 
h 



E, 



(9) 



h 

V04/ 



(36) 



The 5x5 matrix T is presented in appendix B. Using the explicit expressions for 
the elements a, 6, d, a', 6', d' (see appendix B) we find 



det(r-^) 



2c^z^ 



Ca z^, 



(37) 



with 
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C5 



8ci 



^ 



^ 



1 + 6$2 + $^ 



$2)2 



$2 



2\ 2 



$2 



^ /I + $2 



C4 



2i^c2 



^ 
^ 



2\ 2 



$2 



^ 



1 2 



1-^ 



$^)^ $4 



i + *y (i-$2)6 

C5 is negativ and hence the largest eigenvalue in fl36l) is 



£;, 



(9) 



-C5 



C4. 



1 + ^ 
1-^ 



(38) 
(39) 

(40) 



As we need to know the largest eigenvalue of the transfer matrix for all positive 
values of J/ J', we have to also consider the space of odd states. The appropriate ansatz 
for an odd eigenstate is 



0«(7*) = <t>il*i+ <t>2 ill + h Hi 73 74 + </'4 7i7273 

+ 0-1 72 + 0-2 ill + 0-3 ^727374 + 0-4 7l7274 • 



(41) 



Here the products of Grassmann variables in the first four terms have total 
momentum q mod tt, while the products in the last four terms have total momentum 
— g mod TT. As the momentum q is conserved (mod vr) in our model, the eigenvalue 
equation resulting from 



^d\*d\ Ti/2(7*,A) e 
is block diagonal, 

/0±i \ 

0±2 

0±3 

V 0±4 / 



-A*A 



0«(A*) 



.{«) 



Er uii 



(42) 



r± 



/0±i \ 



E, 



(n±) 



(43) 



P±2 

t>±3 

\ 0±4 / 

Explicit expressions for the matrices 7+, 71 are to be found in appendix B. The 



two largest eigenvalues y E, 



.{u+) 



Eq turn out to be identical, 



E, 



iu+) 



E, 



(«-) 



-C4. 



(44) 



As a first application of the results obtained so far in this section, we identify 
the transition line between ordered and disordered states of the AIK. As is shown in 
section VB of reference [5] a transition from a high-temperature disordered phase to an 
ordered state at low temperatures occurs when the even and the odd eigenstates become 
degenerate as the temperature is decreased. Thus we have to compare the eigenvalues 
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E^ = (— C5 + a/cj — C4 y and E^ = C4 as functions of the temperature T in the different 
regions of J/ J'. 

1) J I J' < 1: Using the exphcit expressions for C4 and C5, equations ( l39l) . ( l38l) . one finds 
after some algebra that E^ > E^ for all g 7^ 0. For q = E^ and E"^ become equal if 



l-tp 
l + ^p 



20 



or J = — In 



Tc sinh^(^) 



2 cosh( 



2J' 



(45) 



This defines the transition line between the ferrimagnetic and paramagnetic phase 
of the AIK (see figure H] below) which was already determined by Kano and Naya. 

2) J = J', isotropic kagome model: Here Eg is larger than Eg for all momenta 
q and for all temperatures including T = 0. This reflects the fact that the ground state 
of the Ising model on the kagome lattice is disordered. 



3) J/ J' > 1. Using the explicit expressions for C4 and C5 again, one finds that 
Eq/Eq attains its maximum at g = 0. As T decreases to zero, Eq/Eq -^ for 
g 7^ 0. However, for g = 0, the ratio Eq/Eq — )■ 1 as T — )■ 0, implying that there is 
long range order in the ground state of the AIK. 

These facts will be discussed in more detail below after we have determined the 
correlation function between the middle spins. 
Since Eq > Eq for all g > 0, the free energy per site of the AIK is 



0<q<7r/2 

In figure [2] we show the entropy per site that results from this expression. 



Figure 2. Entropy S per site of the AIK model. J' has been set to unity here. 
Temperatures (from bottom to top): 0.05, 0.3, 0.5, 0.75, 1.0, 1.5, 2.5, 10.0 in units of 
J'. As T — ?► 00, the entropy S reaches its hmit In 2. As T — > 0, 5 decreases to for 
J < J' while a rest entropy of i In 2 remains for J > J' (sec introduction). 
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Obviously, the isotropic point J = 1 plays an exceptional role in the 
thermodynamics of the AIK. 

4. Correlation function between middle spins 

For a system with M rows of spins the correlation function between two middle spins 
at sites m and n {m < n) of the same row is given by 

tr { T,%' al^al | 

X(m, n) = —^ ^ (47) 

tr { T,f } 

where T1/2 is the transfer matrix (128|) defined in section 3. If the state |3) with 
maximum eigenvalue of the transfer matrix T1/2 is non-degenerate this reduces to 

X(m.n) = ^^|g^ (48) 

in the thermodynamic limit M — > 00. Again, the further development in this 
section follows closely the method devised by Schultz, Mattis and Lieb [9] in section 
IV of their work. Applying the Jordan- Wigner transformation, equation (fT2|) to the 
product a^ a^ we have 

al,al = ihlKe'-^^--^^'>\ (49) 

where 

6J = h] + hj and ih] = h] - bj . (50) 

Expanding the exponential on the r. h. s. of equation f l49|) we obtain for the 
correlation function fl48l) 



X{m,n) = (0,|(z6^)(6^+i) i^bl^,) ■ ■ ■ (bU) i^K-i) (6^)|</>,)/(</>,|</>,). (51) 

Here, l^g) is the ground state of — lnTi/2 in the even subspace in the fermionic 
representation (c. f. equation fl33|) . 

7r/2 

!</>.) = n \^9),^ 

where 

- Uq) ^i^\A + 4.^1) + 04(g) ^1,4,^^^,4,10) (52) 

(|0) : fermion vacuum). 

The fermion operators b\^, b\^, b\^ and 64^ act in the intervals [0,7r/2[ , [— 7r/2,0[ , 
[—re, — vr/2[ and [re /2'n[ , respectively. Our choice of the even ground state in equation 
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(IHTj) implies that by evaluating this expression we will obtain the correlation function 
in the disordered region of the phase diagram where \(j)g) is the ground state. Indeed, 
we expect to find the most interesting results in this region. 

Since \(j)g) is the ground state of a system of noninteracting fermions, cf. section 3, 
the expectation value in equation flSTl) can be evaluated by the use of Wick's theorem, 
which requires the knowledge of the contractions of all pairs of operators of the product 
in equation ( ISTl In appendix B the coefficients 0,/(g), z/ = 0, ■ ■ ■ , 4 are found to be real 
for all q. This suffices to show that 







(53) 



((/.,|6J6f|0,) = O and (0,|6pri0,) 
for all pairs of sites j ^ I. Denoting the remaining contractions by Gj^i^i, 

G,,,+i = (0,|z6pf+,|0,)/(0,|0,), (54) 

we can thus cast the result for the correlation function flSTl) into the form 



X(m,n) 



'-^m+l,n— 1 ^m+l,n 



G 



n— l,m+l 



G. 



71— 2, n 



G, 



n— 1,71 



(55) 



In appendix C we find that Gj^i = 0, if j and / are both even, or both odd, and that 
Gj^i = 5j+i,«, if J is odd. Furthermore, the remaining elements Gj^i are found to depend 
only on the difference I — j. Hence, denoting 

G2j,2i+i = a-i-j (56) 

we get for the correlation function between middle spins at the sites 1 and 2n + 1 



X(l,2n+1) 



1 











■ 








Go 





ai 


■ 


fln-l 








1 





■ 








a_i 





flo 


■ 


an-2 











1 


■ 

















■ 1 








Q'-n+l 





a'-n+2 ■ 


■ 


flo 
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ai 



O'Tl-l 



(57) 



O'-n+l CL-n+2 " " " ^0 

Thus, for the AIK, the correlation function between the middle spins is represented 
by a Toeplitz determinant, just like the spin-spin correlation function of the square 
lattice Ising model, cf. reference [9]. 
In appendix C we show that 



where 



dz 1 1 

2m z z^ 



P2+P1Z' 



Poz' 



1-1=1 



\P2+P1Z-^ +P0Z~ 



(5^ 



Po 

Pi 

P2 



(l_$2)^Vl/2 



(1 - ^Y [(1 - *')' (1 + ^^) + 4$' (1 - ^) (1 - *^ 



[(1 



$2)2 ^ 



2$2(1_^)2]^ 



(59) 



This expression for the elements of the Toeplitz determinant (157|) agrees with the 
corresponding expression in reference [18] apart from a sign difference in pi, which is 
probably due to a misprint in [TS] . 

As is described in detail in Chapter XI of reference [16], for large n the Toeplitz 
determinant fl57|) can be evaluated analytically by applying a technique that was 
developed to solve Wiener-Hopf sum equations. We apply this technique here, but 
refrain from presenting any details apart from the remark that, as in the case of the 
square lattice Ising model, we have to deal with the case of index u = —1 here (see 
reference [16], p. 250). The result takes the form 



X{l,2n + 1) = {-)" R{n + l)x. 



(n) 



Asymptotically, 



:i-4)(i-^7 



2Nni/4 



and 






dz 
27ri 



.n-l 



1 — Z Zb 1 — Z ^ Zb 



z z 



1-1=1 




z~^ z 



1 .-!• 



(60) 



(61) 



(62) 
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(Here and in the sequel the relation symbol " ~ " implies that corrections to the 
right hand side of the expressions vanish in the limit n — )■ oo.) The real quantities z^ 
and zb are the roots of the the quadratic form p2z'^ + piz + po (see equation (l58l) ). 



za,b 




(63) 



Since pi < 0, the roots za and zb are positive and zb < I < za for all values of the 
parameters $ and \1/ outside the ordered region. However, depending on the values of 

K = J/T and K' = J'/T either ® zb < z^^ or ® zb > zj^. 



® Zb < z/ 

We evaluate the integral in equation ( l62l) by deforming the integration contour such 
that it encircles the branch cut of the integrand which extends from zb to 2^^, 



(64) 




In the limit of large n we obtain (see equation ([60 



X(l,2n 



[1 - z%){l - ZaZb] 
l-zi" 



1/4 




Tin 



(65) 



® Zj^ < Zb 

Proceeding in the same manner as in case ® we find 

^'^dt 






t 



n-1 



1 — t Zb Zb — t 



-1 IT 



A. b Z A fj ^ A 



and hence in the limit of large n 



X(l,2n + 1) 



:i-zini-z-/) r ZAZB Y 

{1-zbz2')^ \zaZb-1J 



nl/4 



Ann^ 



-B 



(66) 



(67) 



We remark here that for the expressions (!65l) and (1671) to be valid, the distances n 
must be large enough so that the condition 



n:^ n= \\n{zAZB] 



(6^ 



is satisfied. 

The plot of n in figure [3] indicates that in large regions of the phase diagram, 
in particular for J > 1, the correlation function x(l,2n + 1) attains its asymptotic 
behaviour only for rather large n. 
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Figure 3. n {J' — 1); as J — > 0.91, n — > since zazb — > in this limit, cf. equation 

mi- 



5. Results and discussion 



The most remarkable feature of the phase diagram, figure HJ is the disorder line Jd{T), 
which separates the regions ® and ® defined in the previous section from each other, 
see reference [23] and references therein. It is determined by the equation 



$2 = (tanh(J7T))2 = tanh( J/T) = ^ . 



(69) 



Figure 4. Pliase diagram (T in units of J' , J' is set to unity). FM: ferrimagnetically 
ordered region, full line: critical line Jc{T); CD, ® disordered regions, separated by the 
disorder line Jd{T) (dashed line). 
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On this line 

= -. (\z\ = 1) (70) 

Then it follows from Eq. f l58|) that in this case aj = 5j^_i and hence x(l; 2r2 + l) = 0, 
i.e. the middle spins are uncorrelated on the disorder line. 

5.1. Unfrustrated Ising like region 

In region ®, the asymptotic behaviour of the correlation function is given by 

_ 2n 

X(l,2n + 1) ~Cr^-i, A-i = i^^^^, Ci > (cf. equation (ESI)). (71) 
'n 2 



This suggests that in this region the AIK model behaves like a unfrustrated Ising 
magnet in the disordered phase, cf. reference [16], Chapter XL In fact, for J = the 
AIK model, cf. figure [H reduces to the antiferromagnetic Ising model on the square 
lattice whose edges are decorated by the chain spins, which are decoupled from each 
other. Tracing out these spins yields an effective ferromagnetic coupling between the 
middle spins with coupling constant K^ = — ^ lncosh(2i^'). Furthermore, for J = 
equations ( l63l) and ( 159|) yield 



zb = and zl^ = -— = smh^(2Ka). (72) 

Pi 



Then, the expression fl65|) for the correlation function between the middle spin of the AIK 
model reduces to the known expression for the correlation function between spins along 
one of the diagonals of the ensuing square model in the disordered region, cf. reference 
[16], Chapter XI.2. 

To check the quality of the asymptotic expression f l65|) for the correlation function 
x(l, 2n + 1) we have calculated the determinant in equation ( l57l) numerically. In figure 
Elthe asymptotic result ( 16 5 p and the numerical result are compared after the exponential 
factor e~^"/^^ has been divided out from both results. As is expected on account of the 
condition ( l68l) . for small and moderate distances 2n the agreement is worse if J is close 
to the disorder line where fi — )■ oo. 

On the critical line Tc{J), cf. equation ( H5|) . we find za = ^ and 



«. = -/—--■ J4 '"'" • (73) 

^ J 2TTi z zi \ z^ 1 - Z Zb 

\z\=l 

This is very similar to the corresponding expression for the square lattice Ising 
model, equation (5.1) in Chapter XI of reference [16]. Applying the same method as 
McCoy and Wu, we arrive at 

X{2n + 1, 1) ~ [(1 - zbT{1 - zD] ^'/'^ -^ {A : see reference [16]) (74) 



n 
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Figure 5. Comparison of the asymptotic (line) and the numerical results (symbols) 
for x(l, 2n+ 1) for J = 0.7, n = 0.3 (left) and J = 0.825, n = 17.5 (right) after division 

by e-2"/^i, see text (J' = 1). 



where 



zb 



1 + 2$^ - $4V 1 - 2$' 



$ 



4\2 



A = 0.645 ■ ■ • 



(75) 



(l + 8$2 + 10$4 + 8$6 + $8)2 ' 

In summary, the results presented here indicate clearly that in region ® the features of 
the AIK model differ only quantitatively from those of an unfrustrated Ising model with 
antiferromagnetic nearest neighbour interaction. The middle spins of one row of the 
model are correlated ferromagnetically. The appearance of a ferrimagnetically ordered 
region for small coupling J and low temperatures fully agrees with this picture of an 
unfrustrated model. 



5.2. Chain dominated region 

The features of the AIK model change qualitatively when the coupling J is increased 
across the disorderline into the region @. The disordered phase extends down to T = in 
this region, and the correlation function of the middle spins changes from the expression 
in equation ( I7T|) to 



X(l,2n+1) 



with 



Co- 



^2' 



e «2 



\ia{zB) 



(76) 



C2 > (see equation ([6 



This implies that all the middle spins in one row of the AIK model are 
predominantly correlated antiparallel to one another. In Appendix D we present a 
simple perturbative evaluation of the correlation function which is valid for J' <^ J 
and J' /T <^ 1. In this limit, the correlations between the middle spins of one row are 
mediated by the spins of two adjacent chains which are correlated antiferromagnetically 
and to which the middle spins are weakly coupled. This mechanism leads quite naturally 
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Figure 6. Comparison of the asymptotic (line) and the numerical results (symbols) 
for x(l, 2n + 1) for J = 0.91, n = 0.1 (left) and J = 1.5, n = 54.4 (right) after division 
by e-2«/?2 (J/ = i)_ 



to a negative correlation between the any pair of middle spins of the same row. 
As we have mentioned in the Introduction, for J > J' the middle spins decouple from 
the chain spins in the limit T — >■ and become therefore completely uncorrelated. To 
see how this happens as the temperature is decreased to zero we expand za and zb 
(cf. equation ( l63l) ) for low temperatures, i. e. e = 2e"^'^'/^ <^ 1, ^2^ = e"^^'^"'^')/^ -C 1, 

z2' = l- 2^2' + ^^2"' + (^fe-'e, ^2^ ^2"'), 

ZB =1-2^2" 

With these expansions, equation fl67|) yields 
X(l,2n + 1 



k2' + 0{^2'e,^2'e,^2') 



(77) 



2^2^ 



2 g"2n/6 



(78) 
vanishes in 



Although the correlation length ^2 tends to infinity as T — )■ 0, x(l) 2ra + 1) 

this limit since the prefactor ^^^ vanishes. 

To compare the asymptotic result ( 167|) for x(l, 2?t,+ 1) in the parameter region ® with the 

exact result (!57|) we have evaluated the latter numerically as in the other regime. The 

comparison is shown in figure El As in region (D, the agreement between the asymptotic 

and the exact results becomes worse for small and moderate distances n as n increases. 



5.3. The point J' = J = 1,T = 

From our discussion of the rest entropy of the AIK model in the introduction, it follows 
that in the limit T — )■ this model behaves irregularly at the isotropic point J' = J = 1. 
Here, we wish to discuss the behaviour of the correlation function between the middle 
spins as the isotropic point is approached from different directions in the phase diagram. 
We fix a direction by choosing a fixed value for the parameter a in the equation 



J — J' = aT, 



(79) 
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which represents straight hnes with different slopes in the phase diagram that 
intersect at the isotropic point. For instance, a = —In 2 is the slope of the phase 



boundary T = Tc{J 
at this point and a 
parametrising 



J') at J' = J = 1, a 



-^ ln2 is the slope of the disorder line 
is the line perpendicular to the abscissa in J' = J = 1. By 



$ = 1 -e ^ = 1 

we obtain from equations (!59|) . 



zb 



1 2 



e v^ where s 



4 + 6s 



-4a 



(80) 



Zb 



(2 



" 4 + 6s " 


2 


2 

_2-s_ 



This yields for the different cases: 
s = 16: on the phase boundary x as in equation (^^, 



16 > s >4: 


in region ® 


s = 4: 


on the disorder line 


4> s>0: 


in region @ 


s = 1: 


on the perpendicular 


s = 0: 


on the line T = 



X as in equations (165|) . ( ITT]) 



X 







ZB = 1/49 
Za > Zb 
Za = Zb = ^ 
X as in equations f l67|) . fl76|) zb > zj^ 
X as in equation (176|) ^2 — 1-25 

X = cf . equations (!60|l . z^^ = zb = I 
fl6T]) and (l62]) as on the 

disorder line 



6. Summary and outlook 

In this paper we have adapted the transfer matrix technique devised in reference [9] 
for the solution of the square lattice Ising model to the case of the Ising model on the 
anisotropic kagome model (AIK). The isotropic kagome model is fully frustrated, but 
as the coupling J in one of the three lattice directions is varied from zero to infinity, the 
model changes from the unfrustrated square lattice model to a set of antiferromagnetic 
chains. We have calculated the pair correlation function of middle spins, i. e. of spins 
on intermediate rows between the chains with coupling J, (i) asymptotically exact, (ii) 
numerically and for the case J ':$> 1 (iii) by a perturbative method. We find that the 
phase diagram of the AIK model is divided into two sections by a disorder line Jd{T). 
On this line, the correlation between the middle spins vanishes exactly. The following 
statements need to be corroborated by a calculation of the correlations between the 
chain spins, but this was not a subject of the present paper. For J < Jc{T) the model is 
ferrimagnetically ordered: we expect the chain spins and the middle spins to be ordered 
ferromagnetically but antiparallel relative to each other; in the region JdT) < J < Jd{T) 
we expect short range order of the same kind. (Short range ferromagnetic order between 
the middle spins is found in the present paper.) In summary, for J < Jd(T) we expect 
the AIK model to behave like an unfrustrated Ising model. 

At the disorder line, the correlation between the middle spins changes qualitatively. For 
Jd{T) < J there is no long range order at any temperature. The three methods (i - 
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iii) show that the pair correlation function of the middle spins is negative throughout 
this region. It decays exponentially with a correlation length ^2(^) that diverges in the 
limit T — )■ 0. However, as T — > the amplitude of the correlation function becomes 
proportional to ^2(T)~'^, i. e. the middle spins become uncorrected as T — ;■ 0. Since 
the chains are antiferromagnetically ordered at T = 0, the loss of correlation between 
the middle spins was to be expected. 

To complete the picture of the spin correlations of the AIK model various correlation 
functions in addition to the middle spin correlation function need to be calculated. We 
expect to be able to perform the necessary calculations by applying the transfer matrix 
technique devised here. 
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Appendix A. Working with Grassmann variables 

Here we collect rules to be followed in transforming operators from their fermion rep- 
resentation, e. g. expressions (I3l) . . . (E]) for the transfer matrices Tq ...T3, into their 
representation by Grassmann variables and of the rules that were important in our ma- 
nipulations of such operators. We adhere to the terminology of reference [22] . 
If T{b\ b) is a normal ordered fermionic operator {b\ b: fermion creation and annihila- 



tion operators), then T(7*,7) := T 



, where 7 are Grassmann numbers and 7* 



their complex conjugates, is called the "normal form" of T, and T{'-f*,'-f) := T(7*, 7) e"^ '^ 
is called the "matrix form" of T. {b^ , b and 7*, 7 can be vectors of fermion operators 
or Grassmann numbers, respectively.) 
Gaussian integral: For an arbitrary N x N matrix H we have 



Yl [ d-f*J-f^ e~^nHnmlm = det H. (A.l) 

(In the exponent, automatic summation over the indices is implied.) 
Gaussian integral with external sources: for an antisymmetric 2N x 2N matrix R: 



I 



A X* \R M \+[ 7 7* 

dX*dXe V^V 




2 1 7 7* 1-^1 



7 



±Vdeti?-i g V^ / (A.2) 
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Here, 7 and A are A^- dimensional Grassmann vectors. The determinant of an 
antisymmetric matrix is always positiv. The sign of the square root depends on the 
order of the Grassmann integrations. It is irrelevant for our purposes. 

Trace: 

tr \f] = [ d-fd-f* e^*^ T(7*, 7) (A.3) 



Product of fermionic operators: 

A-B = C: C{-f*-f)= [ d\*d\A{-f*X)B{X*-f)e-^'^. (A.4) 



In our case all products in question can be evaluated explicitly since the integrals that 
have to be performed are gaussian (cf. equation (129|) ). 

Appendix B. The matrices T, 7± 

As the matrices T, 7+ and 71 play a central role in our analysis of the AIK model, we 
list their elements explicitely here in terms of our basic variables ^ and $ (cf. equation 
([8])). T is a real symmetric 5x5 matrix and T± are real symmetric 4x4 matrices. 
Defining T and 7± by 

T=V-T, Tpm = V-Tpm (B.l) 

we have 

Til = 1, T12 = a, Tvi = d', TiA = \/2b', Ti5 = ad' + b'^, 

T22 = 0.^ + a'^, 723 = ad' — b"^, T24 = v2{a'b' — ab), 
T25 = a' {old! + 6'2) - a!b'^ + dla^ - 2abb' , 

%z = d^ + d'^ , fsi = V2{bd + b'd'), 
T35 = d! {a'd' + b''^) + a'd^ - d'b^ + 2bdb' , 

%^ = ad + b^ + 2b'\ %5 = V2 (b' {a'd' + 6'2) 

+a'bd - d'ab + b'{ad + 6^) 

Tss = {a'd' + b'^f + a'^d^ + d'^a^ 



-Ab'^ad - 2a'd'b^ + Aa'b'bd - Ab'd'ab 



(B.2) 



and 



7±ii = a, 7±i2 = —b, 7±i3 = ad' =f bb' , T±u = Tab' — ba, 
7±22 = d, 7±23 = ^db' — bd' , 7±24 = a'd ± bb' , 
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7^33 = ad"^ T '2bb'd'^ + db'\ T±m = Tah'd! ± da'h' - b {a'd' + b'^) , 
7^44 = ab'^ + da'^ ± 2bab'. (B.3) 

In terms of the parameters $ and \1/ the quantities v, a, d, b and a', d', b' are given by: 



2\ rIS4 



16 (1 - c,'^') (1 + ^Y - 4 s/$'^ (1 + $^) + 4 (^^ - c,2) $ 

(1 -$2)^(1 -^2)2 



with 



a = 2Cg {Cq + 1)$2(1 - $4)(1 - ^)2(l _ ^2)/^g^ 
(i = —a\ 



Iq^TT+q 



b = a 



Cq + l 



den 



1 + $4 _ 4$^^ + ^2 ^ 0*^2)2 ^ W(^^ _ <|,^^)^(<|>^ _ ^ 



,2\2 



One sees immediately that ad — b'^ = 0. This identity has been used in estabhshing 
the matrices in equations (jB.2p and (IB.Sp . 



2,T,\2 



^f((l + $')(! + ^2) _4$2^N 
4Cq(l - <l>=^^)(^ - $2)($2(^^ _ ^^,2 ^ ^^^(^^ _ ^4-^A 



d' = -a'\ 



q-^q+n 



b' 



den 



$ 



2n2 



1 + *^)^(1 + $^) + 2<I>2((1 - ^2)2 _ 2* - 2*'^ 



(B.4) 



With this information exphcit expressions for the eigenvalues and the eigenvectors of 
the transfer matrices that are needed in the main part of the paper can be obtained. We 
refrain from presenting the rather lengthy explicit expressions for the coefficients 0i . . . 04 
of the eigenvector ( l33l) but we note certain symmetry properties of these coefficients that 
will be used to simplify expressions to be derived in the next appendix. First we observe 
that 



and furthermore we find the following symmetries of i 






H^ 



^2^ 



03 -^ 



(B.5) 
>!, 03 and 03 as functions of g: 



(B.6) 
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Appendix C. Calculation of the elements of the Toeplitz determinant 

In this appendix we present the main steps of the calculations that lead to the result 
fl58|l for the elements aj of the Toeplitz determinant fl57|) in the main text. We need to 
determine the expectation value (cf. equation ( l54l) ) 

Gm,n = -j^{'PgWJ)l\'Pg) ) ^here Ng = {(pg\(pg) 

= TTi(t>g\ - bm K. - {hm bn)^ + bl^bn + (C bn^ \(t>g)(,'^ " ^m.n) 

+ (^-^(0,|2 6L&„|0,)-l^ 5^m,n, (C.l) 

cf. equation (150|) . Let us first consider the expectation value 

fm,n = -rjr{4>g\bmbn\(pg). (C.2) 

With the representation 



N 

0<q<TT/2 

of the annihilation operator at site m by the momentum space operators big ... 64^ 
introduced in the main text (see paragraph below equation ( !52l) ) we find after a straight- 
forward evaluation of expectation values of the type q{<t>g\b^q b,^q\(()g)g, (/i, z/ = 1, . . . , 4) 

U,n~ Yl ;^ [- 2 sin(g[m-n])(l-(-r+") 00(01 + 02) 

0<q<n/2 ^'^ 

- 4 cos(g[m - n]) ((-)- - {-)-) ^0o] . (C.4) 

Here 

K, = g{<Pg\<Pg)q = 20o + <Pl + <Pl + <Pl (C.5) 

(The momentum label q has been suppressed on the r.h.s.) In equations (1C.4P and (JC.SP 
use has been made of the identity (IB. 51) . Note that fm.n is real and hence 

J^{4>g\bi bU(pg) = J^{4>g\bm bn\^g) = fm,n ■ 

9m,n = Tr((f^9\bln bn\<Pg) (C.6) 



Similarly, 



1 v^ 1 



0<q<TT/2 ^^ 



2 



J2 ir 2cos(g[m-n]) <j 0^ + 0^ + '^3 



2 



+(-r+"(0^+02 + |) 

- 2sin(g[m - n]) ^ [(-)- - (-)"] (0i - 0^) | . (C.7) 
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Obviously gm.^n is real too, and with equation (IC.Sp it follows immediately that 
gm,m = 1/2; SO that the last term in the expression fIC.ll) vanishes. On account of the 
symmetries ( ]B.6p . the sums over q in Eqs. ( 1C.4|) . ( JC.71) can be extended to the full 



Brillouin zone — vr < g < vr, 

N ^L^ AN ^^ 



0<q<TT/2 -Tr<q<n 



Then, 



[i-(-r'i7n-^ + 7.[(-r-(-)i7;-^ 



'^mjU ^ \ Jm,n ' gm,n)\^ Om,n) 
= \[^-{-r-"']ln-m + \ 

with 

^' ^ iV 5Z JT A^'^^^ '^^^ ^^'^^ ~ '^^^ ^°^(^^) ~ 200 sin(g/)] 

and 

ii = I E }^ 4| ['^o ^<<^) + ^(<^i - <^2) sin(g/)]. (C.8) 

— 7r<g<7r ^^ * 



-m 



For the elements of the determinant (I55p this yields 

G2m,2m' = , G2m+l,2m'+l = 

G2m-l,2m' = l2{m'-m)+l — 72(m'-m)+l — ^m'- 

G2m,2m'+1 = 72{m'-m)+l + 72(m'-m)+l — ^tm'-m • (C-9) 

Further simplifications of these expressions are possible on account of the following 
properties of the coefficients 0i, 02 and 03 that were found by analysing the explicit 
expressions for these coefficients with computer algebra: They have a common 
denominator Do{q) which is invariant under the operations q ^ —q and 5 — ;■ vr + g, 

M<i) = ^,, y = 1,2,3. (c.io) 

Only two of the four coefficients -Do? -Di, -D2 and D3 are independent. The relations 
between them are 

- -^cMq) - ^s,[D^{q) + D^iq)] = D,{q) 

and 

-^s,DM + \c,[Dr{q) + DM] = \{DM - D^iq)). (C.ll) 

Using these relations one finds 

a; = 5,,o (C.12) 
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and, introducing p = 2q a.s the momentum variable and denoting D^{^) = D,^, 



«^=/v E 



— 7r<p<7r 



Dl-\{D,~D,f 



cos{p[j + 1]) 



+ Do pi - D,) sm{p[j + 1])} j^ —2. (C.13) 

In the continuum hmit, A^ — )■ oo, and with z = e*^ as the integration variable, this 
leads to the expression ( 158|) for aj. 

Appendix D. Perturbation expansion in J' /T 

As we have mentioned in the Introduction, for J > J' the chain spins of the AIK model 
are ideally ordered antiferrromagnetically in the ground state, and as a consequence 
the middle spins are completely decoupled from the chain spins. For weak interchain 
coupling J' <^ J the order of the spins along the chains will prevail over a long distance 
for low but finite temperatures T <^ J, and the coupling between the chains will remain 
small as long as J' ^ T. In the sequel we will calculate the correlation function between 
the middle spins of row / of the kagome lattice, x(2?^ + 1, 2n + 1) = ( cr2m+i,i o'2n+i,i ), 
perturbatively with respect to J'. (In a^y the second subscript / denotes the row / of the 
lattice.) The hamiltonian of the AIK model couples the middle spins of row I to chain 
spins of the rows / — 1 and / + 1 , 

J 0'2m+l,l ■ [o'2m+l,l-l + 0"2m,i-l + <^2m+l,l+l + 0"2m,i+l] • (D.l) 

In second order in J' /T one obtains for x'- 

X(2m + l,2n + 1) 

J' ( J' 

= 2 ( — )^([0"2m+l,i+l + a2m,l] [o"2n+l,Z+l + O-2n,l])j'=0 + O ( ( — )^ 



,J', 



2 Xid-M„g(2n - 2m) + Xid-ising(2n - 2m + 1) 



Xid-i.ingW = ^-j" (tanh(-) ) . (D.3) 



+ Xid-isi„g(2n - 2m - 1)J + O ( (^)M (D.2) 

Here Xid-isingl'^) is the correlation function of the one- dimensional Ising chain, 

V ? 

Hence 

X(2m + 1, 2n + 1) = - A_^ (^tanh(-) j + O (^{jfj . (D.4) 

In the limit considered here, the correlation between any pair of middle spins of 
the same row is negative for all distances. It is mediated by the interaction with the 
neighbouring spin chains. 
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